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Abstract 

The observed optical/UV spectra of most Seyfert galaxies are much 
redder than expected from a stationary accretion disk. Two explanations 
of this efect are probable: (i) the observed spectrum is strongly contami- 
nated by starlight (ii) the accretion disk is not stationary in these objects. 

The standard accretion discs are known to be thermally and viscously 
unstable over a certain range of temperatures. In the inner disc regions 
there may develop radiation pressure driven instability, which is possibly 
related to the rapid variability detected in AGNs in the UV range. In the 
outer disc develops the ionization instability, similar to that in the cata- 
clysmic variables, but operating on much longer timescales. Due to this 
process the spectrum of the accretion disc differs from that of a stationary 
one. 

We examine the accretion disc vertical structure model in order to 
determine the range of the ionization instability. We derive the radial de- 
pendence of the accretion rate in a non- stationary disc between outbursts, 
and we calculate the effective temperature profile and the corresponding 
spectral slope. 

The predicted shape of the spectrum is different from a stationary 
case but it does not seem to reproduce the observed spectra of Seyfert 
galaxies. We therefore suggest that the starlight contribution is most 
probably responsible for the spectral shape in red Seyfert galaxies and its 
contribution extends to shorter wavelengths than usually adopted. 



1 Introduction 

Accretion disc in Active Galactic Nuclei are subject to thermal - viscous 
instabilities which operate in the inner disc regions. The radiation pres- 
sure driven instability (Lightman & Eardley 1974) is developed in the 
models, which assume the viscous torque proportional to the total pres- 
sure. It operates in the innermost disc part, while further away, where 
the disc is gas pressure dominated, there exists a partial ionization zone. 
Therefore these regions are unstable due to the same mechanism as in 
cataclysmic variables (e.g. Smak 1984, Lasota et al. 1995), but in much 
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longer time scales. Finally, the outer parts of the disc may be gravitation- 
ally unstable (Hure 1998). 

Here we consider only the ionization instability. Its radial extension de- 
pends on several parameters, like the mass of a central black hole, viscosity 
parameter a and external accorotion rate Mext- Using the paramctrization 
given in Siemiginowska et al. (1996) we estimate for the inner and outer 
radii of the unstable zone that Ri„ ~ 2hQRschw and Rout ~ l^QQRschw 
for M = IO^Mq, a = 0.03 and Mext = O.lMEdd- 

At any given radius the local thermal equilibrium may be represented 
by the chracteristic shape of the S curve on the Te// — S (or M — S) 
plane. The upper and lower branches of the positive slope describe the 
stable configurations and the middle part of the S curve describes the 
unstable solution. When the external accretion rate corresponds to the 
unstable branch the disk undergoes oscillations between the upper hot and 
the lower cold branch. Faint phase corresponds to the slow accumulation 
of the material and increase of the disk surface density; the accretion rate 
and surface density corresponding to the upper end of the low branch are 
characteristic for this phase. 

2 Vertical structure model 

The basic equations that describe the disc vertical structure are the equa- 
tion of viscous energy dissipation, the hydrostatic equilibrium, and equa- 
tion of energy transfer. The last equation takes into account the pres- 
ence of convection which carries non-negligible fraction of energy. The 
frequency-averaged opacity k (Rosseland mean) includes the electron scat- 
tering as well as the required free-free and bound-free transitions. The ta- 
bles are from Alexander, .Johnson & Rypma (1983) for log T < 3.8, from 
Seaton ct al. (1994) for log T > 4.0 and the value of opacity is interpo- 
lated between these two tables for intermediate values of the temperature, 
as in Rozanska et al. (1999). 

The temperature dependence of the Rosseland mean opacity changes 
for different temperature ranges: for T < 10^ K the power law function 
is: K ~ pT~^, for T > 10"' K the power law index is equal to -3.5, while 
at the intermediate temperature range, corresponding to partial hydrogen 
ionization, the trend is opposite and power law index is about 5-10 (see 
Kato et al. 1999). The opacity changes affect the cooling rate in the disc, 
which is important in determining the equilibrium solutions. Moreover, 
the presence of ions causes a large reduction of radiative transport 
efficiency and therefore excited convection results in higher opacity. 

3 Results 

In Figure 1. we plot local stability curves for several different radii. The 
critical turning points on each curve, indicated as ai, 02 and 03, represent 
the maximum surface densities achieved during the faint phase of the 
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FIGURE 1. The theoretical stabihty curves for M = 10 M©, a = 0.03 and three values 
of radius r: r = IQORschw (solid line), r = 350Rschw (dashed line), and r = WOORschw 
(long dashed line). The accretion rate is given in Eddington units. 

ionization instability cycle. The critical mass transfer rate is Mcritif) = 
M(E™,,(r)). 

In Figure 2. we show the radial dependence of accretion rate in the 
faint phase, derived from our vetical structure model. The accretion rate 
is in Eddington units and the radius is in Schwarzschild radii. This de- 
pendence is well fitted with the power law function: m(r) ~ r^'^. 

The effective temperature of the accretion disc scales with radius as 
T(r) ~ (M(r)/r^) and for the stationary accretion disc this relation gives 
T(r) ~ r~°-''^, while in our model we obtain T(r) ~ r~°-"^ . The ap- 
proximate spectral slope is given by the formula ~ u^'^^^, where p 
is the index in the radial dependence of the temperature (Kato et al. 
1999). In the first case it results in Fi, ~ v^^^. In the latter case it 
becomes Fi, ~ however as the outer disc temperature is not very 

much smaller than the inner disc temperature, the resulting spectrum 
looks practically like a single black body. None of these, however, is ob- 
served in the Seyfert galaxies, for which the optical/UV spectral indices 
are in the range from -1.9 to -1.2 (e.g. Edelson & Malkan 1986). For NGC 
7469 the spectral slope in the optical band is ~ —1.69 (see also Nandra 
et al. 1998), while for the faint state of NGC 4151 it is ~ -1.7 (Lyutyi & 
Doroshkievich, 2000). 

Therefore our conclusion is, that even if the instability processes may 
affect the emitted optical spectral slope, nevertheless the non-negligible 
effect of starlight component should be taken into account. 

In Figure 3. we show the example of the starlight spectrum, taken to 
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FIGURE 2. The radial dependence of accretion rate (vertical structure model parameters: 
M = IO^Mq, a — 0.03) . The accretion rate is given in Eddington units. The dashed line 
represents the power law function m(r) ~ r^'"^. 



be the spectrum of the nuclear bulge of M31 galaxy, as in Kuraszkiewicz et 
al. (1997). It is compared with the power law spectrum of the index -1.9 
( e.g. NGC 4051, NGC 3516, Mkn 231). We conclude, that the starlight 
contribution may be responsible for this spectral slope, however it must 
extend to shorter wavelengths. Clearly more young, blue stars are needed 
in this case to account for this shape of the spectrum. 
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FIGURE 3. The starlight spectrum (solid line) compared with the power law vFy ~ v 
(dashed line). 
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